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Abstract 

Inspired by an old idea of von Neumann, we seek a pair of commuting operators X, P which are, 
in a specific sense, "close" to the canonical non-commuting position and momentum operators, x,p. 
The construction of such operators is related to the problem of finding complete sets of orthonormal 
phase space localized states, a problem severely constrained by the Balian-Low theorem. Here 
these constraints are avoided by restricting attention to situations in which the density matrix 
is reasonably decohered (i.e., spread out in phase space). Commuting position and momentum 
operators are argued to be of use in discussions of emergent classicality from quantum mechanics. In 
particular, they may be used to give a discussion of the relationship between exact and approximate 
decoherence in the decoherent histories approach to quantum theory. 

PACS numbers: 03.65.-w, 03.65.Yz, 03.65.Ta 
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I. INTRODUCTION 



A. Preamble 

At the heart of quantum mechanics is the canonical commutation relation between posi- 
tion and momentum operators, 



Physically, this relation corresponds to the fact that measurements of position and mo- 
mentum depend on the order of the measurements. It is this relation, in essence, that is 
responsible for the key differences between classical and quantum mechanics, since in clas- 
sical mechanics, measurements of position and momentum can be made in such a way that 
the order makes no difference. It follows that any account of the emergence of classical be- 
haviour from quantum theory must reconcile these two very different aspects of classical and 
quantum theory. A typical point of view is that classical mechanics emerges only at a very 
coarse-grained level and for sufficiently coarse-grained samplings of position and momentum, 
their non-commutativity makes little difference Q|. 

Still, one wonders whether there is a deeper or more precise way of reconciling the non- 
commuting quantum operators with their commuting classical counterparts. Indeed, this 
clearly troubled the founders of quantum theory, since von Neumann addressed the issue in 
his 1932 book. Mathematical Foundations of Quantum Mechanics P]. He noted that when 
we make observations of a macroscopic system, we are in fact able to make observations of 
position and momentum simultaneously (although imprecisely, of course). This suggested 
to him that the measurements we make do not in fact correspond directly to the usual 
operators x and p, but to some other operators, X and P, say, which commute 



and which must in some sense be "close" to the original operators x, p. Such a pair of 
operators could be particularly useful in bringing a degree of precision to discussions of 
emergent classicality. The aim of this paper is to discuss the construction and utility of such 
operators. 



[x,p] = ih 



(1) 



[X,P]=0 



(2) 
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B. Classical Imprecisions 



To what degree do macroscropic measurements fix the position and momentum opera- 
tors in quantum mechanics? When we make macroscopic measurements of, say, a particle, 
there will be imprecisions Ax, Ap in the specifications of position and momentum. These 
imprecisions may be "small" to classical eyes but they will typically be very large compared 
to the quantum scale. It will be useful for what follows to get a quantitative idea of this. 
Suppose that the position is measured to a precision of 10~^m and the velocity to within 
10~^m/s (both extremely precise specifications, by macroscopic standards). For a mass of, 
say, 10~^kg, we then have. 



This means that when we specify the phase space location of a classical system in a way that 
is very precise to classical eyes, at the quantum scale it could be in any one of about 10^^ 
phase space cells. There is therefore a considerable amount of freedom at the quantum level 
to redefine the position and momentum operators without making any noticeable difference 
to macroscopic observations. Perhaps within this freedom there is the possibility to find 
position and momentum operators which commute. 

C. Motivations 

What are the motivations for constructing commuting operators which are close to the 
position and momentum operators? An important indicator of emergent classicalitvis ap- 
proximate diagonality of the density operator p (often referred to as decoherence) [3]. For 
example, for a point particle interacting with a thermal environment, it may be shown that 
after a short time scale, the density operator approaches the form 



where f{p,q) is a non-negative function and {ippq) are a set of phase space localized states 
(such as generalized coherent states) This indicates that it is approximately diagonal 
in both position or momentum. (It cannot be exactly diagonal in both, since they do 
not commute.) This means that there is negligible interference between different values 
of position and momenta and, loosely speaking, they may then be treated as if they are 




(4) 
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classical. The statement is imprecise, however, since there is still some interference, so the 
variables are only imprecisely defined. 

One would like to be able to make a more exact statement about diagonality of p and this 
is possible in terms of the commuting operators X, P. For example, using the eigenstates 
\ipnm) of X, P one could construct a density operator of the form 

P ~ ^ fnrn\'4'nm) {lpnm\ (5) 
nm 

which is exactly diagonal in X and P, indicating there is exactly zero interference between 
different values of these quantities. One can then choose the coefficients fnm to make p' as 
close as possible to p. The question of the degree of approximate diagonality is therefore 
shifted to the question of the closeness of the density operator to a pseudo-classical density 
operator p'. 

A more cornprehensive approach to emergent classicality is the decoherent histories ap- 
proach llll |. There, the central object of interest is the decoherence 
functional, 

= Tr {PaM ■ ■■PaAtl)pPa[ih) ■ " ■P<(tn)) (6) 

The histories are characterized by the initial state p and by the strings of projection operators 
Pa{t) (in the Heisenberg picture) at times ti to t„ (and a denotes the string of alternatives 
ai ■ ■ - an)- Intuitively, the decoherence functional is a measure of the interference between 
pairs of histories a, a'. When 

ReL'(a,aO=0 (7) 

for a 7^ a', we say that the histories are consistent and probabilities p{a) = D{a,a) obeying 
the usual probability sum rules may be assigned to them. (Typically, the physical mecha- 
nisms producing consistency actually cause the stronger condition -D(a, a') = for a 7^ a' 
to be satisfied, which is referred to as decoherence.) One can then ask whether these prob- 
abilities are strongly peaked about trajectories obeying classical equations of motion, and if 
they are, we can say that the system is emergently classical. 

For histories in which the projections Pa^i^k) are onto positions at different times of 
a point particle interacting with a thermal environment, the decoherence functional, like 
the density operator, is approximately diagonal 0, Q]. The approximation is typically 
exceptionally good, but still only approximate. Again one wonders whether more exact 
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statements can be made. Indeed, it has been conjectured that approximately consistent 



121. There are a number 



histories can be in some sense distorted into exactly consistent ones 
of ways in which a set of histories could be distorted: one can change the initial state, the 
times of the projections, the widths of the projections or the operators who spectrum is 
projected onto. The possible existence of the commuting variables P, X suggests a particular 
way of distorting the histories so as to make them exactly decoherent. 

Each position at time t is, in the Heisenberg picture, a function of the canonical pair, 
x,p, so Xt = ft{x,p). The projections onto Xt at different times do not commute, since 
7^ in general. The decoherence functional is therefore not diagonal in general, but 
can be approximately diagonal if the system is coupled to a thermal environment. 

Now suppose we consider projections onto the variables Xt = ft{^, P) at different times. 
Under reasonable dynamics, Xt will be close to Xt as long as x,p are close to X,P. The 
operators Xt do commute at different times so all the projections commute and, as is easy 
to see, the decoherence functional will be exactly diagonal. So exact decoherence is achieved 
by shifting the operators x,p to the commuting pair X,P. Furthermore, it is known that 
the probabilities for histories of positions Xt are typically peaked about classical evolution. 
The probabilities for the commuting variables Xt will therefore have the same property if 
the Xt are close to Xt- The question of approximate decoherence and emergent classicality 
is therefore shifted to the question of the closeness of the old and new operators. 

These then, at least in outline, are the reasons why a commuting pair of position and 
momentum-like operators may be useful for discussing emergent classicality. We turn now 
to the construction of such operators. 



D. Von Neumann's Construction 



Von Neumann outlined a prescription whereby the commuting operators X and P may 
be constructed [^. This involved first taking a discrete subset \m,n) of the coherent states 
\p,q), with one state per cell of size 27ih (a von Neumann lattice). He alleged that these 



states are complete (this was later proved 



13 



He then stated that they may be 



orthogonalized using the Schmidt process to produce an orthonormal set \ipnm)- From these. 
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he constructed position and momentum-like operators 

X = '^na\ijjn-m){^nm\, P = m'^^\tjjnm) {^nm\ (8) 

nm nm 

(where a is a constant with the dimension of length). These operators clearly commute. He 
argued that these new operators are indeed "close" to the old ones, in the sense that 

(^„^|(X - X) Vnm) ii^nmliP " P) Vnm) < K^fl" (9) 

Von Neumann's calculations imply that the constant K is about 1, 800 (but he thought more 
detailed calculations could give a smaller value). 

However, von Neumann's prescription is at best at sketch of how this works and he 
certainly did not give full details (such as the explicit form of the {ipum))- Furthermore, we 
now know a lot more about phase space localized states than was known in 1932, and, as will 
be described below, there are obstructions to constructing such states. These obstructions 
do not necessarily apply to what von Neumann did, but nevertheless, it is still interesting 
to revisit his ideas from a more modern perspective. 



E. A General Approach and the Balian-Low Theorem 

We start with a more general statement of the problem. We consider a set of states of 
the form 

l^nm) = Unm\i^) (lO) 

where {ip) is a fiducial state and Unm is the unitary shift operator, 

Unm = exp (^^no,P - ^mbx^ (11) 

A particularly interesting case is that of a von Neumann lattice, in which case b = 2TTh/a 
and the translations in the p and x directions then commute, and we have 

Unm = (-1)"" exp ('-nap] exp ( -t'^x] (12) 



h J \ a ^ 
There is then one state per cell of size 27ih, as in von Neumann's case. It is of interest to 
find a fiducial state in Eq. lfTIH) such that the states are complete and orthonormal, that 
is, 

'^\^nm){^nm\ = 1 (13) 
nm 

{i>nm\iJn'm') = ^nn'^mm' (14) 



Fiducial states leading to states satisfying these properties are easily found and we will 
exhibit a set below. There is, however, a crucial difficulty. According to the theorem of 



Balian and Low 
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17, 



18|, if the three properties (jl(J|) . (|13|) and (jl4j) are satisfied, then 



the fiducial state has either (Ax)^ or (Ap)^ infinite, so is not phase space localized. If 
we used such states to construct commuting operators X, P as von Neumann did, then at 
least one of the averages in Eq.Q would diverge, and there would be no sense in which 
the new operators are close to the old. (Note that von Neumann claims to have used the 
Schmidt procedure to construct his orthonormal set, which one would not expect to produce 
states satisfying Eq. (|Tn|l . so his construction does not necessarily fall foul of the Balian-Low 
theorem) . 

The problem of constructing orthonormal phase space localized states is one of great 
interest in a number of fields so some effort has been expended in finding ways around the 
Balian-Low theorem. Zak has proved some interesting results in this area. In Ref [lol, he 
showed that the coherent states restricted to a von Neumann lattice |m, n) obey a sort of 
orthogonality relation if the usual inner product (m, n\m\ n') is averaged over a single phase 
space cell. It is not yet clear if this result can be used to produce commuting position and 
momentum operators. He has also considered complete orthonormal sets of states which 
are localized in position, but double-peaked in momentum (so localized in p^, but not in p. 
j'rom these one can construct commuting operators X and P^, which are "close" to x and 
2ol | . This is tantalizing close to the goal of this paper, but not quite there (and also suggests 
that the p —p transformation plays a crucial role in the Balian-Low theorem). Many of 
these and similar results are proved using the so-called kq representation, a technique which 



is particularly well-adapted to these problems 

To avoid the Balian-Low theorem, one has to drop one of the three requirements (|TU]). (fT^ 
and (fT^ in order to get phase space localization. For the purposes of this paper, which is to 
construct useful commuting position and momentum operators, the orthogonality Eq. (fT^ is 
essential. We will therefore explore the possibility of dropping the other two requirements. 
Eq. (|Tn|l . the requirement that the states be obtained by translation of a single fiducial state is 
mainly for practical convenience, so there is no harm in relaxing this as long as the resulting 
states are not unmanageable. 

More significantly, we will relax the requirement of completeness, Eq. ()13|) . The motive 
behind this is as follows. We are interested in using commuting operators to discuss emergent 
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classicality. In practice, this means that we are only concerned with the physical situation 
in which the density matrix of the system has undergone a degree of decoherence. This 
means that it is approximately diagonal in both position and momentum, or equivalently, 
its Wigner function is reasonably spread out in phase space. The density matrix is therefore 
insensitive to the fine structure of phase space and it seems reasonable to suppose that the 
physics could be well-described by a less than complete set of states, if they are carefully 
chosen. Of course, this is a quantitative matter that needs to be checked in detail and we 
will do this. 

For the purposes of constructing commuting position and momentum operators, we will 
actually use a construction slightly more general than the one indicated by von Neumann 
in Eq.(jH)). In particular, we will look for commuting operators X, P of the form 

X = Y^^nEnm (15) 

nm 

P = Y^P^E^^ (16) 

nm 

Here, X„ and Pm are c-numbers and Enm are projection operators localized onto a region of 
phase space labelled by n, m which will consist of more than one 27r^-sized cell. They are 
exclusive 

EnmE-n'm' E^m^nn'^mm' (!'''') 

and exhaustive 

Y,Enm = l (18) 

nm 

We will also insist that they are obtained from unitary shifts of a single projector E, 

Enm UnmE (Unm) (l'^) 

Here, Unm shifts from one cell to the next, so is of the form Eq. (fTT|) . with a, b chosen so that 
the translation in the position and momentum directions commute, but ab > 2TTh. These 
three conditions are the natural generalization for projection operators of the requirements 
firU]) . (fT^ and (mi), and the original case is obtained with the choice 

Enm ~ |V'nm)(V'nm| (^^) 

The quantity TiE is a measure of the number of 27r/i-sized cells projected onto, so clearly 
TtE = 1 in the pure state case, but Tr£' 1 in our case, as we will see. 
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One might have thought that this more general construction could avoid the Balian-Low 
theorem, since the restrictions on Enm are in fact weaker than the restrictions (fTUI) . (fT^ 
and (|14p for pure states. This is not in fact the case. We will prove a modest extension 
of the Balian-Low theorem which shows that there is no phase space localized projector E 
satisfying the three requirements Eqs.(|17|). (|18|) and (|19|). However, what we will do is find 
an "almost" complete set of phase space localized pure states from which we can construct a 
projector E that satisfies the exhaustivity condition Eq. ()18|) to a good approximation when 
acting on sufficiently decohered density operators. From this we can construct commuting 
operators X and P may with useful properties. 



F. Earlier Work 



Finally, we briefly mention a related approach. An earlier attempt to construct commut- 
ing position and momentum operators was considered in R,ef . [2^ . This construction involved 
doubling the orig inal Hilbert space and then using operators defined on this enlarged space. 
(See also Ref. 23[ for similar ideas') . The resulting theory is essentially the same as 't Hooft's 



deterministic quantum theory 



24| . However, this is no longer standard quantum theory. In 



the present work, by contrast, we stay within the framework of standard quantum theory. 



G. This Paper 

In Section 2, we briefiy summarize some known properties of the Wigner function which 
help to make precise the idea that the state is sufficiently spread out in phase space. We also 
briefly note that the Wigner function naturally suggests an alternative method of deflning 
commuting position and momentum operators. In Section 3, we prove a modest general- 
ization of the Balian-Low theorem for projector operators. In Section 4 we introduce an 
orthornormal set of phase space localized states that are "almost" complete. We then use 
them to construct a set of phase space localized projection operators Enm which are almost 
exhaustive. In Section 5 it is shown that the incompleteness does not matter if the density 
operator of the system is reasonably spread out in phase space. In Section 6, we use the 
projection operators Enm to construct commuting position and momentum operators and 
compute the "distance" between these operators and the usual canonical pair, as in Eq.Q. 
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In Section 7, we compare the probabilities for X and P with those for the usual position and 
momentum operators and find them to be close. We summarize and conclude in Section 8. 



II. SOME PROPERTIES OF THE WIGNER FUNCTION 

It will be useful for the rest of the paper to briefly summarize here some aspects of 
quantum mechanics in phase space and the Wigner function. Most of this is standard 
material and may be skipped by the informed reader, except the brief comments at the end 
of this section. 

The Wigner representation for a density operator p (or indeed a wide class of operators) 
is defined by 

Wip, l) = ^Jd^ piq + l^,q- ^0 (21) 

with inverse 

p{x, y) = j dp e>(^-^) W{p, ^) (22) 

Many calculations involving operators are usefully expressed in the Wigner representation. 
For example, 

Tr {AB) = 2Txh j dpdq Wa{p, q)WB{p, q) (23) 

where Wa{Pi q) and Wb{Pi q) are the Wigner functions of A and B. 

We will be interested in later sections in the behaviour of the density operator or Wigner 
function in a simple model of the decoherence process. We take the simplest case of a single 
particle coupled to a thermal environment in the limit of high temperature and negligible 
dissipation, with no external potential. The master equation for the density matrix p(x, y) 
is, 

dp ih (d'^p d'^p\ D 2 



at 2m\dx^ dy^J ^2^"^ P ^^'^^ 
where D = 2m'jkT. In the Wigner representation, the corresponding Wigner function obeys 
the equation, 

dW__pdW ^d'^W 
dt m dx dp^ 

Evolution according to the master equation Eq.(j2D) tends to produce approximate diago- 
nality in position and momentum. In the Wigner representation, this appears as a spreading 
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out phase space. Indeed, using Eq. (j25j) it may be shown that 

(Ap)^ = 2Dt+iAp)l (26) 



where 



(j{x,p) = -{xp + px) — {x){p) (28) 



evaluated in the initial state. In particular, for long times, the phase space spreading behaves 
according to 

(^p)'(^^)' . (JUL) e (29) 

h \ h J 

This means that any initial state becomes spread out in phase space on the (typically very 



short) timescale {h/'jkTY^'^. (See Ref. 25| for similar calculations). 

In the case above, a free particle, the spreading continues indefinitely. However, for a 
bound system with dissipation, equilibrium is eventually reached. For a simple harmonic 
oscillator at thermal equilibrium, for example, the ratio of thermal to quantum fluctuations 
is 

iAp){Ax) ^ kT 

At typical laboratory temperatures and frequencies that are not unrealistically fast (for 
macroscopic systems), this number can be very large, of order 10^° say. These estimates will 
be relevant later on. In brief, they show that the density operator becomes very spread out 
in phase space (and hence slowly varying) very readily. 

This much is known material and will be used later. However, we now note that the exis- 
tence of the Wigner representation suggests an alternative method for constructing commut- 
ing position and momentum operators. It is well-known that when the Wigner function is 
sufficiently spread out, it becomes positive and may be regarded as a probability distribution 
for the variables p and q (which clearly commute, since they are numbers, not operators). 
The variables p and q do not correspond exactly to the operators p and x, but they are close 
when the Wigner function is spread out. In fact, it is easy to see that we have the following 
correspondences: 

qW{p, g) ^ ^ (xp + px) (31) 
pW{p,q) ^ ^{PP + PP) (32) 

11 



That is, multiplication by g or p in the Wigner representation corresponds to the operation 
of anticommutation with x or p on the density operator (and it is easy to show that these 
two operations on p commute). The point here is that these operators on p are operations 
on the space of density operators which have no counterpart in terms of operations on pure 
states. This is therefore not a route to producing the desired pair of commuting operators 
envisage by von Neumann. 



III. AN EXTENSION OF THE BALIAN-LOW THEOREM 

We first show that there is no projection operator satisfying the three properties Eqs. ljlSp . 
(jl9|) and (j'iUj) . This is an almost trivial extension of the proof of the Balian-Low theorem 



given by Battle 



id- 



We consider the object Tr(£'xp), which we assume exists. (If it does not, i.e., is infinite, 
this means that E has infinite dispersion in either x or p, so is not phase space localized). 
We have, from the three properties of -Emn, Eqs. ffT7|) -(fT ^ . 

Ti{Exp) = ^ Tr (ExEmnP) 

mn 

= J^T^r (ExUmnEULp) 

mn 

= J2T^^iEUmn{x + na)E{p + bm)Ul^) 

mn 

= Tr {E^m~n{x + na)E{p + hm)) 

mn 

= "^Tt {E_rn,-nXEp) 
mn 

= Tr (xEp) 

= Tr {Epx) (33) 

Or in other words, 

Tr(E[x,p]) = (34) 

which, via the commutation relations, implies that TrE = 0. Since E > 0, this means that 
E = 0. So no projector satisfying the three properties exists, if one insists that they be 
satisfied exactly. Hence, as indicated we will relax the conditions in what follows. 

We note in passing that this simple result has implications for the Balian-Low theorem in 
the pure state case. It may seem that one could avoid the Balian-Low theorem by dropping 
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the requirement Eq.(fTU|) and requiring that the states iV'nm) are generated from more than 
one fiducial state. However, one could then use those fiducial states (assuming they are 
orthogonal) to construct a projection operator E satisfying the properties Eqs. (fT7j) - (fT^ and 
the above result shows that the Balian-Low theorem is not in fact avoided. 



IV. AN ALMOST COMPLETE SET OF ORTHONORMAL PHASE SPACE LO- 
CALIZED STATES 

We now show how to construct an almost complete orthonormal set of phase space local- 
ized states, from which we can construct the projector E. Our starting point is the set of 



states considered by Low 



iJnmi^) =a-'2h{x-na) e^™/'^ (35) 



where h{x) a window function on [—^a,^a]. These clearly satisfy the three conditions, 
Eqs. fll()|l . ()13j) and ()14j) . The Fourier transform of these wave function is 

from which it is easy to see that (Ap)^ diverges because ipnmip) goes to zero like 1/p for 
large p which is not fast enough. Differently put, the window function in Eq. ljH^j) causes 
the derivative of the wave function to involve the 5-functions 6{x — na± |a), which are not 
square-integrable. These properties are fully in line with the Balian-Low theorem. 



However, inspired by a suggestion of Zak 



3,i 



it is easy to see from Eq. (jHH|) that we may 



take linear combinations of these states which fall off like and therefore have finite 
(Ap)^. In particular, the states 

i^nli^) = {'^n,2m{x) + V^n,2m+1 (x) ) (37) 

have this property. One can also see this in configuration space: they vanish a.t x = na± |a, 
and the offending 5-function appearing in their derivative therefore causes no problems. 
They are not complete, consisting of just "half" the states in the momentum direction, and 
indeed the states missed out are the states 

Xnlix) = {'lpn,2m{x) - iJn,2m+l{x)) (38) 
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and these have infinite dispersion in p. 

But it is not hard to see that the "remainder" states Xnm decay hke 1/p in momen- 
tum space, so a further process of "halving" is possible to produce more states with finite 
dispersion. That is, we define 




with new remainder states 

X^nl = (xi'L(^) + xl'L+i(^)) (40) 

So the set of states V'im, V'im, Xnm is complete and orthonormal, but the states ^jJnm have 
infinite dispersion. 

We can continue in this way to define the sequence of states 

'^nm' ^ ^ (^Xnjm+l ~ xiijrnj (41) 

and 

(K+l) _ ±_ ( (K) (K) \ . 

for K = 1,2,3, ■■■ . If we truncate the sequence at some finite value of K, K = N, say, then 
the set of states ipnm, 'tpnm ■ ■ ■ i^nm together with the remainder states Xnm are orthonormal 
and complete. We therefore have the completeness relation 

N 

E E ^i'SM^niiyr + E xi'2i^)xi'2iyr = s{x - y) (43) 

K=l n,m n,m 

In some sense, "most" of the states, namely the -^im', have finite dispersion, and "some" 
of them, namely the Xnm have infinite dispersion. In this way, as N increases, the infinite 
dispersion anticipated from the Balian-Low theorem is pushed into a progressively smaller 
set of states. (An interesting question is whether the limit N ^ oo may be taken in any 
meaningful or useful way, but we will not pursue that here). 

Since each state is a linear combination of the ipnm, one may also derive the following 
general formula, 

2^-1 

i^fn^ = 2-^/^ E # ^n,2^m^j (44) 
j=0 
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where 



and 



"""-{-ly if 0<j<2^-^-l (45) 



cf^-{-iy if 2^-^ < j <2^-l (46) 



We also have 



This may also be written, 



(47) 



2^-1 

Xi^^ = 2-^/^ J] (-1)^- V'„,2-.+, (48) 

j=0 



(K) ^ 2-^/2 ii ! U (49) 



The states \'4^l!^) are not all obtained from a single fiducial state, since we have 

\^l^) = Ui^\r) (50) 

where 

Ui'S = t/„,2-^ (51) 

and 

IV'^) = Ko) (52) 

There are therefore fiducial states for the set iprm, K = 1, ■ ■ ■ N . 

At some length, one can compute the averages and dispersions of p and x in the fiducial 
states {i/j^}. One obtains 

= (53) 



a \ 2^ 
{x)k = ' (54) 

(Ax)^ = g (56) 

(Note that {p)k is not zero. Because there is more than one fiducial state, it does not 
appear to be possible to shift the states so that = in the fiducial states without 

spoiling orthogonality.) 
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The construction described above is concisely summarized as follows: For each set of 2^ 
lattice points in the momentum direction, there are 2^ — 1 states with finite dispersion and 
just 1 state (the remainder state), with infinite dispersion. 

Consider the following simple example to illustrate the construction. Consider the 8 
lattice points m = 0, 1 • ■ • 7 in the momentum direction, so A'^ = 3. Then the 8 states which 
depend only on these points are 





1 

i'ipnfl + V'n,l) 




<! = 


{tpn,2 + 1pn,2) 




Yn,2 — 


^ (V'n,4 + V'n.s) 




^b^^^ - 
Yn,3 — 




(57) 



V'S = ^ ii^nA - 1pn,5 - 1pn,6 + Ipnj) (58) 



^n,0 = ^ (^n,0 - V'n,l + ^n,2 " ^n,3 " ^n,4 + ^n,5 " ll'nfi + V'nj) (59) 
Xn,0 = ^ (V'n,0 - 1pn,l + 1pn,2 - 'tpn,3 + tpnA - 1pn,5 + 1pn,6 " Ipnj) (60) 

There are 7 finite dispersion states and one infinite dispersion state. It is easy to sec that 
they are orthonormal and form a complete set on these 8 lattice points. Note also that 
widths of the states increases with K, but this is in some sense offset by the fact that 
the states become progressively sparser. 

So far we have been working with a complete set of states. We now come to the specific 
form of the proposal to relax the requirement of using a complete set of states (in the 
construction of commuting operators P and X): wc quite simply drop the remainder states 
Xnrn, which have infinite dispersion, and use only the incomplete set , K = 1,2 - ■■ N, 
which have finite dispersion. That is, in each set of 2^ lattice points in the momentum 
direction, we use 2^ — 1 out of the 2^ states. 

It seems likely that this approximation will be valid for sufficiently large N, for suitable 
density matrices. We will show in the next section that reasonable results can be obtained 
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using the ifj^ alone. Here, we briefly look at the remainder states Xnm to see under what 
conditions they may be dropped. Note first that the states ipi^{x) have the property 
that they vanish at the end-points of the intervals, x = na ± a/2 (as they must, so that 
their derivative does not have a 5-function). Not surprisingly therefore, the remainder 
states Xnm{x) become narrower and progressively more concentrated about the end-points 
as — > oo, as one can see from Eq.(|49p (since they are in some sense making up for the fact 
that the i/j^{x) vanish at the end-points, and the whole set of states is complete). This 
suggests that, under suitably coarse grained conditions, the behaviour of theses states at 
single points will become insignificant. We will see this more explicitly below. 

We may now use the almost complete set of states to construct the desired projection 
operators Enm localized on phase space cells. We choose the phase space cells to have 2^ 
lattice points in the p-direction and 1 lattice point in the x-direction. We have found 2^ — 1 
states with finite dispersion in each of those cells. We therefore take the projector E to be 

N 2^-^-1 

^ = E E (61) 

K=l m=0 

It is a sum over all states depending only on the lattice points to 2^ — 1 (in the momentum 
direction). Importantly, this is an exact projection operator which localizes onto a region of 
phase space - it satisfies 

E^ = E (62) 

exactly. The projection operator for any other cell is easily obtained by unitary displacement, 
using steps of size 2^ lattice points in the momentum direction (and single steps in the x- 
direction) : 

Enm = UiZ^E{Ui^jy (63) 

These projectors clearly satisfy the exclusivity condition, Eq. (jl7j) . exactly, but do not satisfy 
the exhaustivity condition Eq. ()18|) since we have 

+ E \xl!2){xi!2\ = i (64) 

nm n,m 

But, as we have argued, we expect the x terms to be negligible under suitable conditions so 
the projectors Enm should be almost exhaustive 



J^Enm^l (65) 



nm 
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The approximate nature of this property may not in fact matter for many practical apph- 
cations. The phase space projector onto a large cell F in phase space is defined by 

n,m€r 

The projector onto the region outside F is then defined to be 

Er^l-Er (67) 

so they are trivially exhaustive, E-p + E-p = 1. The key point here is that the remainder 
states dropped in the construction of E have infinite dispersion so they do not naturally 
belong in the construction of a phase space projector for a finite region of phase space. 
We will need some further properties of E. We have 

TiE = 2^-1 (68) 

which means it projects onto a phase space region of size (2^ — l)(27r/i). The object E/TrE 
may be thought of as a density operator and we can compute averages and variances to see 
the properties of E. We have 

- ^ (69) 
^a^(.--i) (70) 

and 



(AP)I-E^ pi) 



2^ 

K=l 

2^+i7r2/i- 



2 



3a2 

where we show only the leading terms of large N. We also have 



(72) 



{x)e = (73) 



(Ax)| = g (74) 

Since (p)^ 7^ 0, it is in fact useful to perform a simple translation in momentum and define 
a related projector E' with all the same properties as E except that Tr(p£") = 0. We will 
use this in what follows. 
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The construction of an exact projector with the above properties is the main achievement 
of this section and will be used to construct the commuting X and P operators below. 

Finally, we make two minor remarks. First, note that position and momentum enter in 
the construction in very different ways. However, the constant a is arbitrary, so may be 
tuned to make the width of the projector E arbitrarily small in either the x or p direction. 
Also, the states Eq. ()35p and ()36|) are a Fourier transform pair, so we could easily interchange 
them and start with a set of states that have perfect localization in p, instead of x. It would 
of course be of interest to find a construction in which x and p entered on an equal footing. 

Second, we note for comparison that Omnes has made extensive use of phase quasi- 
projectors of the form 



where \p,q) are phase space localized states, such as the coherent states j9|. These have 
proved very useful for discussing emergent classicality in quantum theory. However, in 
contrast to the projectors constructed here, these are not exact projectors, since they obey 
the approximate relation 



(and so Pr and 1 — Pr are only approximately exclusive). It would be of interest to revisit 
some of Omnes results using the exact projectors constructed here. 

V. VALIDITY OF APPROXIMATE COMPLETENESS 

Now we come to a crucial check of our approach, which is to determine the conditions 
under which working with an approximately complete set of states gives reasonable results. 
The density matrix p of the system satisfies Trp = 1. If the states \ipi!^r!) for K = 1, - ■ ■ N 
are approximately complete, then we should have 




(75) 




(76) 



N 




(77) 



K=l nm 



which is the same as 




(78) 



nm 
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We check this. It is most useful to exploit the Wigner representation, Eqs.(j2I]),(I221) together 
with the property Eq. lj^ . so we have 

EE<^-H^-) = 2^^EE dpdqW^^\p,q)W,{p+ m,q + na) (79) 

K=l nm K=l nm 

where Wp is the Wigner function of p and W'^^'^ is the Wigner function of Now the 

crucial step is to approximate the discrete sum over ra, m with an integral over continuous 
variables, p = 2^^^TTh'm/a, q = na and we have 

ETT. / 2^^^7ih . f dpdq , _ , 

Wp{p + m,q + na) ^ / ^K+i^f^ Wp{p + p,q + q) (80) 

nm 

^ (81) 



This approximation is valid as long as the Wigner function of p is slowly varying over phase 
space volumes of size 2^{2t[K). Since K runs up to A^, we require slow variation on a scale of 
size 2^{2ttK). This is typically the case for density matrices that are sufficiently decohered, 
as we saw in Section 2. We now have 

N N ^ 

EE^^i^^H^i™^) - (82) 

K=l nm K=l 

= 1-^ (83) 

Hence we do indeed get a result close to 1 as long as is sufficiently large. 

One other related result is worth recording here since it will be used in the next section. 
Any density operator satisfies the relation 

/ ^U\q,k)pU{q,k) = l (84) 

where [/(g, k) is the unitary shift operator in phase space. This is easily proved using the 
Wigner representation above, since we have 

{x\U\q, k)pU{q, k)\y) = j dp e^^^^-^^ Wp{p - k, ^ - q) (85) 

and integrating the right-hand side over k and q yields 27rM(a; — y). 

One would expect that, for sufficiently slowly varying Wigner functions, a discrete version 
of the result Eq.(jHH) would hold. So suppose instead of the operator U{q,k), we use the 
operators Un^ defined in Eq.d^H). We then have 

{x\{Ui^^ypUil^\y) = jdp e^^^-y^ Wp{p - ^^m, ^ - na) (86) 
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Using the same approximation and same steps as in Eq.()8U|) to do the sum over n, m, it 
follows that 

E ^Ul^^ypU^^ ^ i (87) 



2K 

mm 

and therefore 

E E ^UL'SypUL^ - 1 - ^ (88) 

K=l nm 

VI. CONSTRUCTION OF COMMUTING POSITION AND MOMENTUM OP- 
ERATORS 

We now use the projection operator E'^^ to construct commuting position and momentum 
operators, as outlined in the introduction. They are, 

X = E (89) 

nm 

P = J2Pm E'^^ (90) 

nm 

where X„ = na and Pm = m x 2^ {2Txh/a). Clearly 

[X,P] = (91) 



as required. 

We need to determine whether these operators are close to the original canonical pair, 
x^p. To do this we need some measure of distance, || P — p ||. We will define this by 

II P-p f = Tr (^(P-p)V) (92) 

(and similarly for X) where p is the density operator of the system, which we will assume 
is reasonably spread out in phase space. A more general measure of distance would include 
some sort of optimization over p, but, for the class of reasonably decohered density operators, 
the result depends only weakly on p. 
We have 

II P - p f = Tr [{P^ - 2Pp + p2)pj + Tr (^P[p, p]) (93) 

where the operators P are moved to the left in the first term in preparation for inserting the 
explicit form for P below (since [P,p] 7^ 0). It is useful to introduce the notation tr(- ■ ■ ) to 
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denote a trace over the incomplete set of states We may then write 

II P - p f = tr {{P^ - 2PP + p')p^ + c/J + Tr [p[p, p]) (94) 

where 

dl = Tr(pV) - tr(pV) (95) 

In the previous section we saw that Trp and trp are both very close (to order 2"^) and it is 
reasonable to expect the same if p is replaced with p^p, so will negligible and we drop it. 

We also expect the term Tr ^P[p, p]j to be small, since a decohered p will be approxi- 
mately diagonal in momentum. Inserting the explicit expression for P, it is easily shown 
that 

Tr (P[p, p]) = ^-Tr ([p, {Ui'^^)' pUi'^]) (96) 

nm 

Since we are working in the approximation in which p is sufficiently slowly varying that the 
sum over n becomes an integral. It is then easily seen (using Eq.()841 for example) that, in 
this approximation, the object 

^ ^j {^nm)'^ P^nm (97) 

n 

is diagonal in momentum. Therefore the term Eq. ()96|l is zero to the approximation we are 
using. 

Consider the one remaining term in Eg. ()94|1 . Inserting the explicit for for P, we have 

\\P-P\? = J2tr {EUPm-p)'p) (98) 

nm 

= J2'r{E'f{U^Z^ypU^Z^) (99) 



From Eq.(jHZj), we have 



and so 



J2(Ul^ypUnT - ^ (100) 



P-pf ^ {Ap)l (101) 



3a2 

where only the leading term for large is given. 



(102) 
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A similar (and simpler) calculation shows that 



X-x\Y ^ {AxYe (103) 

(104) 



a 
12 



Putting all these results together we obtain 

il P-p II . II X-x II ^ Ch (105) 

where 

C = 2^/2^ (106) 
3v 2 

This result is valid in the approximation of large N and for slowly varying density operators. 

Eqs. flHHj) and ()105|) show that X needs to be chosen to be "large" to get approximate 
completeness, yet "small" for the commuting operators to be close to the original canonical 
pair. It seems likely, however, that there is a range of intermediate values that will meet both 
of these requirements. For example, take X = 20. Then C ~ 10^, safely within the limit 
estimated in Eq.Q, so the difference between the old and new operators will be completely 
invisible to macroscropic observations. The error due to approximate completeness is of order 
, which is about 10~^. So there appears to be a large regime in which the approach works 
well. 



VII. PROBABILITIES FOR POSITION AND MOMENTUM 

The results Eq. (j83|) and Eq. (|l(J5p are only indicative, and the true test of closeness of the 
old and new operators is a comparison of the probabilities for x,p and X, P. 

Consider the probability that the variable X lies in the range Ax, where Ax is the 
interval [nia,n2a\, for some integers ni,n2. The probability is 

p{Ax) = TTipE^^) (107) 

where the projector E^^ is defined by 

(and the loose notation n G Ax means n G [^i,''^2])- The probability for lying outside the 
region Ax is defined using the projector 1 — E^xy so that the probabilities add to 1 exactly, 
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and the approximate completeness discussed earlier poses no problems. The probability for 
P is similarly defined, in terms of a projector E^p defined by 

Both of these probabilities therefore involve the probability Pnm associated with our basic 
phase space cell of size 2^{27rh), given by 

Pnm = Tr{E'^mP) = {E' {U^l^y pU^Z^) (110) 

It is more usefully written in the Wigner representation, 

Pmn = 27ih J dpdq WE'ip, q) Wp{p + Pn,,q + X^) (111) 

where, recall, X„ = na and Pm — rnx 2^{27rh/a). 
The probabihty for the variable X is 

piAx) = Yl J^P-rn (112) 

Since the Wigner function of p is assumed slowly varying, the sum over m may be approxi- 
mated by an integral over P^ regarded as a continuous variable. The p in p + P^ is absorbed 
into the integration and we obtain, 

p{Ax) ~ ^ E / ^^^^ ^) P^^l + ^n,q + Xn) (113) 

Now the p integral may be carried out with the result 

p(Ax) ^aJ2 fdq p{q + ^n, ? + ^n) (114) 

Since E' is phase space localized, one can see that the first part of the integrand is a smearing 
function peaked about ^ = and with width {Aq^ ^ a^/12. It is normahzed to 1 (for large 
N) when integrated over q since Tr£" — 2^ — 1. If we assume that 

» ^ (115) 

then the presence of the smearing function makes no difi^erence and we obtain 

p(Ax) ~ a p{na,na) (116) 

neAx 
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If, as we assume, the density operator is sufficiently slowly varying for the discrete sum to 
become an integral, we obtain 

p{Ax) ^ [ dX p{X,X) (117) 

J Ax 

the usual probability for position. 

Similarly for P, with analogous approximations, we obtain 

^ /27r;i\ ^ f {p\E'\p) 

p{Ap) = 2 (^1 2^ dp p{p + Pm,P + Pm) 

^ ^Nf^\ J2 p(P^,Pj (118) 

where p{p,p') is the Fourier transform of the density matrix p{x, y). Again, for slowly varying 
density operators the discrete sum becomes an integral and we obtain 

p(Ap) ^ / dP p{P,P) (119) 



p 



which coincides with the usual probability for p. 

We therefore find that the probabilities for X, P coincide with those for x, p as long as 
the following conditions hold: 

(i) is sufficiently large that the errors 1/2^ are tolerably small 

(ii) The density operator p is slowly varying on scales of size 2^{2ttK). 

(iii) The widths of the projections satisfy A^Ap ^ 2^(27r/i). 

The key restriction is (ii), the restriction on the density operator. Eqs.()29|) and ()30|) 
indicate that the density operator can easily become sufficiently broad for (ii) to be satisfied. 
For example, for a time-evolving state, (ii) will be satisfied for 

->(^)"^- 

This can easily be extremely short, even for large N. Similarly, for a state close to thermal 
equilibrium, (ii) is satisfied for 

^ » 2^ (121) 

which is again easily satisfied. 
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VIII. SUMMARY AND DISCUSSION 



We have constructed a pair of commuting operators X, P which, at sufficiently coarse 
grained scales, are close (in a variety of ways) to the canonical position and momentum 
operators x,p. These commuting operators offer a new way of defining the relationship 
between approximate and exact decoherence. 

There are two ways in which this programme could be developed and improved. First of 
all, this paper has concentrated on the technicalities of constructing X and P. It would be 
useful to develop more details of the conceptual framework in which they are used to discuss 
emergent classicality. 

Secondly, the present approach works at the coarse grained scales of about 2^ phase 
space cells. Although arguably small to classical eyes, there are some ways in which this 
scale is quite large, and there are indications that some version of the present approach 
should work at finer scales. For example, it is known that the Wigner function (and also 
the P-function) become positive when coarse grained over just one or two phase space cells 
(rather than 2^ cells, as here), a feature that is often taken as indicator of approximate 
decoherence and emergent classicality [22,] • This suggests that there might be another way 
to construct commuting position and momentum operators which does not require such a 
large amount of phase space coarse graining. It would, for example, be of interest to see if 
von Neumann's original suggestion (involving an explicit orthogonalization of the coherent 
states) can actually be made to work. 

These and related questions will be pursued in future publications. 
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